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ABSTRACT 

We present an exactly solvable quantum field theory which allows rear- 
rangement collisions. We solve the model in the relevant sectors and demon- 
strate the orthonormality and completeness of the solutions, and construct 
the S-matrix. In the light of the exact solutions constructed, we discuss 
various issues and assumptions in quantum scattering theory, including the 
isometry of the Moller wave matrix, the normalization and completeness 
of asymptotic states, and the non-orthogonality of basis states. We show 
that these common assertions do not obtain in this model. We suggest a 
general formalism for scattering theory which overcomes these, and other, 
shortcomings and limitations of the existing formalisms in the literature. 



1. Introduction 

Quantum scattering has been an important subject of study since the early days of 
quantum physics. Unfortunately, while we have a reasonable understanding and intuition 
for simple scattering problems, such as single channel scattering, we cannot say the same 
for more general scattering problems such as multi-channel scattering, rearrangement col- 
lisions, field theoretic scattering, problems where bound states appear, and the like. There 
have been many attempts to generalize scattering theory to deal with more complicated 
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cases. However, the literature in this field, though vast, is highly implicit and not con- 
structive. Most authors that have dealt with the problem have carried over the intuition 
developed from the study of single channel potential scattering. This intuition, while quite 
adequate for simple problems, is ill-equipped to deal with more complicated scattering 
problems. Therefore, it is important to examine the common claim by some authors, for 
example Haag [1,2], that their formalism is general enough to encompass complicated scat- 
tering problems, as well as field theory. Unfortunately, most such formalisms are based 
largely on previous results from potential scattering. Furthermore, even when these prob- 
lems are addressed in quantum mechanical scattering, field theoretic scattering remains 
problematical. Many papers, such as the paper by Gell-Mann and Goldberger [3], treat 
field theoretic scattering as somewhat of an afterthought, without much development from 
first principles, or such as the papers by Van Hove [4], treat it for discussion. The 

first clear development of field theoretic scattering from first principles was the seminal 
paper by Lehmann, Symanzik, and Zimmerman [5]. However the LSZ formalism is not 
applicable in many cases, for example, collisions in which stable bound states appear. This 
is, in fact, pointed out by the authors themselves. 

All this leads to the question: how many of our results and assumptions, and how much 
of our intuition can we carry over from simple single-channel potential scattering to more 
complicated scattering situations? To attempt to answer this question, we will construct an 
exactly solvable sector for a quantum field theory. This model has a three particle sector, 
and allows rearrangement collisions. We will use the solutions of this model, along with 
previous results, to point out where the existing formalism has defects and shortcomings. 
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Our model, which we shall call the Rearrangement Model, is an elaboration of the Lee 
Model [6], and the Cascade Model [7], but with extra particles and couplings chosen in 
such a way as to allow rearrangement collisions. The couplings of the model are B <->■ Ccf> 
and D <->■ CO. This model has a sector, which we shall call the Rearrangement Sector, 
B6 <->■ C6cf> <-> D<f>, in which rearrangement collisions can take place. The model can 
be applied directly to physical problems involving rearrangement collisions. We shall, 
however, leave the applications to subsequent work. 

This model is interesting because it a very simple one, and yet contains the essence of 
many phenomena that can take place in an interacting system. It displays the following 
characteristics: 

1. New states can appear, which have no corresponding states in the original 
Hamiltonian. 

2. The thresholds and continuous spectra shift, and the spectra of H and Ho are 
not the same. Furthermore, the continuous spectra are shifted by different 
amounts. 

3. Genuine rearrangement collisions can take place. Yet we have the sub-additivity 
of the spectra: the spectra in the higher sectors is the sum of all the spectra in 
the lower sectors, with possibly additional terms. 

We will construct the solutions of this model, and then, in the light of the solutions we 
have constructed, will examine various assumptions and assertions made in the literature 
about quantum scattering theory. In particular, we will focus on four key points, that of 
the isometry [1] of the Moller wave matrix [8], the normalization [1] and completeness [9] 
of the asymptotic states, and the non-orthogonality of the physical B8 } C8<f), and Dcf> 



4 

states [10,11]. We shall show that these assertions do not obtain in this model. We also 
comment upon the use of the (renormalized) free Hamiltonian in the literature [1,3,9,10,12], 
rather than the correct prescription, which is to use the isospectral comparison Hamiltonian 
(see section 8). 

The plan of this work is as follows. We start with a review of scattering theory in both 
the single channel and multiple channel cases. In the next three sections, we introduce the 
Hamiltonian of the Rearrangement Model, show how explicit solutions can be found for 
this model in the Rearrangement Sector, and verify that the solutions obtained are, in fact, 
solutions to our model. We then show that the solutions obtained are orthonormal and 
complete, write down the Moller matrix and the comparison Hamiltonian, and show that 
the comparison Hamiltonian is isospectral with the full Hamiltonian, but not with the free 
Hamiltonian. Then we calculate the S-matrix of the system in this sector, demonstrate 
its unitarity, and calculate its eigenphases. In section 12, we present a general formalism 
for scattering theory which overcomes the shortcomings and limitations of the existing 
formalisms in the literature. In section 13 we discuss scattering theory and its relation 
to the solutions that we constructed, and to previous work. Finally, in section 14, we 
summarize our work and present our conclusions. 
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2. The Single Channel Formalism 

For our purposes, it makes no difference whether we use the time dependent or time 
independent formalisms of scattering theory. We are concerned with the assumptions and 
results that obtain from the formalisms, and they remain essentially the same in both 
cases. We shall therefore restrict ourselves to the time dependent formalism in the next 
two sections. We follow the treatment of Newton [13] for both sections. 

The discussion in the next two sections is supposed to be very general. In fact, even 
though the method described deals with the non-relativistic region, "the formalism set 
up is such that, provided there exists a consistent relativistic quantum field theory, the 
transition to the relativistic domain is relatively simple" [14]. However, we find that even 
in such a simple model as our Rearrangement Model, these anticipations are not fulfilled. 
This formalism leads to wrong and contradictory results, as will be discussed in section 13. 

We wish to solve the Schrodinger equation 



We assume that this split can be carried out: we shall consider the case of rearrangement 
collisions later. We define four Green's functions 



(2.1) 



We split H into a free Hamiltonian and an interaction Hamiltonian, 



H = H + H'. 



(2.2) 




(2.3a) 



(2.36) 



with the initial conditions 

G + (t) = g + (t) = 0, t<0, (2.4a) 

G-(t) = g-(t) = 0, t>0. (2.46) 

G + and Q + are therefore the advanced Green's functions, and G~ and Q~ are the retarded 
Green's functions. 

These may be solved formally yielding 

G + (t) = -ie- iHot 8(t), (2.5a) 

G-(t) = ie- iHot 0(-t), (2.56) 

g+(t) = -ie- iHt 8(t), (2.5c) 

g-(t) = ie- im 0(-t). (2.5d) 

Let ^o(t) be a state vector satisfying the free Schrodinger equation. The operator G + 
can then be used to express the state vector ^>o(t') for any time t' later than t, in terms of 
its value at t' = t, 

(<') = *G + {t' - t) (t) . (2.6) 
^o(t') then satisfies the free Schrodinger equation for t' > t, and ^>o(t') — > ^o(t) when 

Therefore, 

lim G + (t) = lim g+(t) = -il, (2.7a) 

lim G~ (t) = lim g~ (t) = il. (2.76) 

*— s-o- t->-o- 

Similarly, for t' > t we can write 

*(t')=ig+(t'-t)*(t), (2.8) 
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and for t' < t we have 



(*') = -iG~ (t'-t) *o(<)> 



(2.9a) 



(2.96) 



We now wish to define "in" and "out" states. We start by defining 



(t) = iG+ (t - t') * (*') , 



(2.10) 



whose time development for t > t! is governed by the free Hamiltonian, but which at time 
to was equal to \I/(to)- We now let the time, t', approach ±oo. Then, for the case of 
t — > +oo, we have the "out" state, and for the case of t — > — oo, we have the "in" state. In 
terms of the "in" and "out" states, the equations for ^(t) are 



Note that these are retarded and advanced Green's functions for the whole system. These 
are not the same functions as those that appear in a (time ordered) Dyson series which are, 
instead, time ordered particle propagators. Note also that for every state in the continuous 
spectrum of Hq } and only for such states, these formulae define a corresponding state in 
the spectrum of H . 

If we insert Eq. (2.9b) in Eq. (2.11a), we find 




(2.11a) 



(2.116) 



(2.12) 



where 




(2.13) 



is called the wave operator or the Moller matrix. We can similarly define fi( ). 
Because H' is hermitian, Eq. (2.13) gives us the relation 

M/ in (t) = 0(+) t M/(t), (2.14a) 

and similarly 

^ out (t) = fi(-^(t). (2.146) 

Then, Eq. (2.12) and Eqs. (2.14) give us the relations 

^ in (t) = Q (+ ^Q (+ ^ in (t) , (2.15a) 

*out(<) = 0H t 0(-)^ out (t). (2.156) 

We now consider the possibility that the free states \&i n (t) and ^ \it(t) span the entire 
Hilbert space, i.e. they are complete. From this assumption, we conclude that fi( + ) and 
are isometric, 

0(+) t n(+) = 0(")t (-) = l (2 .i6) 

This does not, however, mean that the fi's are unitary: we cannot conclude from Eq. (2.15a) 
that Eq. (2.16) holds with its factors reversed. 

Furthermore, granted the assumption that the \I/i n and \& ut each form complete sets, 
we conclude that 

fffiW = fi^iio. (2.17) 
When there are bound states in the spectrum of H , we proceed as follows. Let ^o(E } a) 
be the eigenstates of the free Hamiltonian with eigenvalue E, and a be the set of variables 
necessary to remove any degeneracy. Then, the completeness of these states can be written 
as a resolution of the identity, 

poo 

1 = V/ d£^ (£,a)*o T (E,a) . (2.18) 



We again emphasize that we do not know whether the states of Hq are complete, a priori. 
We are simply proceeding under that assumption. We then insert this into the product 
fifit to get 

fifit = / dE V (E, a) ^ (E, a) fit 

poo 

= dE y2v(E,a)W(E,a) 
Jo 

= 1-A. (2.19) 

A is called the unitary deficiency of fi. From the completeness of the set of all states, 
bound and scattering, of H , 

A = E<"X"/ (2-20) 

n 

Thus, A projects onto the space spanned by the bound states of H . If H has no bound 
states, then fiW and fi(~) are unitary. Both H and Hq are hermitian; therefore, the 
hermitian conjugate of Eq. (2.17) gives 

# ftt = ^H. (2.21) 

We now let both sides of Eq. (2.21) act on ty(E,a) to get 

Hotfv(E,a) = E&V(E,a), (2.22) 

which shows that if E is in the spectrum of H but not in the spectrum of Hq then 

fiW(£,a) = 0, 

and so 

fitA = 0. (2.23) 

Thus, the range of the operators fi^ is not the entire Hilbert space. Instead, these 
operators map the whole space onto the subspace spanned by the continuum eigenstates 
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of H . We cannot reach the subspace spanned by the bound states of H , and therefore, 
cannot construct an inverse operator for the whole space. The closest that we can come 
is to use the operators fiW^ which are inverses of fiW on the subspace of states spanned 
by the scattering states of H , and which annihilate the subspace of bound states of H . 

Assuming that the asymptotic states are complete, we construct the S-matrix in the 
following manner. We use Eqs. (2.15b) and (2.12) to write the "out" state in terms of the 
"in" state, 

*out(<) = fiHQ(+)^ in (t), (2.24) 
which defines for us the S-matrix 

S = Q ( -^Q (+ \ (2.25) 

The S-matrix can be shown to be unitary and isometric. See Newton [13] for details; note, 
however, that S is only unitary when O is unitary. (This point is not clearly stated in 
Newton, or in the literature.) 

Some mathematical questions about convergences arise in the above. Conventionally, 
in the Schrodinger picture (the one in which we are currently working), one demands that 
(see Newton [13] for details) 



lim (t) - V>in (<)] 0, (2.26a) 

— > — oo 

lim (t) - Vwt (<)] 0, (2.266) 

-^+oo 

lim Vin (0) => ^ (0) = fi(+Vin (0) , (2.26c) 



where denotes the strong limit. 

We will find, in our model, that if we construct the asymptotic "in" and "out" states 
correctly, these limits will be satisfied; however, the states will not be orthonormal or 
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complete. On the other hand, if we make the usual assumptions of scattering theory, 
namely that the asymptotic states are orthonormal and complete, then these limits will 
not be satisfied. 



3. The Multiple Channel Formalism 

The above formalism is only adequate for simple single-channel cases. For more general 
scattering problems, such as rearrangement collisions, we must generalize the formalism. 
We shall again follow Newton [13]. 

We want to split up the Hamiltonian into two pieces: one piece, H a} that is left when 
the two initial fragments are taken far apart, and the remaining piece, H' a . We can then go 
through the same development of ^(t) from \&i n (t) as above. However, there is a difficulty 
that occurs for the development for the distant future. If rearrangements or break-ups can 
occur, then it is possible that the "channel" Hamiltonian in the future is different than the 
"channel" Hamiltonian in the past. 

The various possibilities for an n-particle system are handled by defining a partition 
of them into k clusters, denoted by a^. Given a partition a, we define H a by allowing all 
distances between clusters to independently tend to infinity. Therefore, H a will contain 
only interactions that are internal to clusters, but none between them. Then, H' a is defined 
by the requirement that H = H a + H' a , and therefore, for any two partitions a and 6, we 
have 

H = H a + H' a =H b + H' b . (3.1) 
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To each partition, there will correspond Green's functions given by 

i^-H a ^G ± a (t) = 16(t), (3.2) 

with the same boundary conditions as Eqs. (2.4). If H a} after removing the kinetic energy 
of the center of mass motion and of the centers of mass of its clusters, has at least one 
bound state, it is called an arrangement channel. When this condition on H a does not 
hold, the channel is not of interest as an initial or final scattering state. If H a has more 
than one bound state, then each of them defines a separate channel, and therefore, in 
each channel the clusters are in a specific bound state but moving freely relative to one 
another. The channel consisting of the entire n-cluster partition is the channel because 
then H a = Hq. 

Now consider the space of each arrangement channel a, which we shall denote by 7i a - 
Then, if a has m fragments, each state in 7i a will have m groups of bound particles. This 
means that unless the channel a is the entire n-fragment arrangement channel, 7i a will not 
be the whole Hilbert space: the ionized eigenstates of H a will be missing. Furthermore, 
as each 7i a is defined by different channel Hamiltonians, H a} the 7i a 7 s are generally not 
orthogonal to each other. In fact, "the complete set of basis functions is not linearly 
independent and, of course, not orthonormal" [15]. 

It will be convenient to define the orthogonal projections P a onto the channel spaces, 
7i a - In other words, we define 

P a = Pa-, Pa = Pa-, PaT^-a = (3-3) 

with the null space of P a defined as the space spanned by the ionized eigenstates of H a . 
P a projects states from the full Hilbert space, 7Y, to the channel spaces, 7i a - Obviously, 
for the n-cluster arrangement channel, we have Pq = 1. 
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We now wish to define "in" and "out" states. We first define an a state, which is a 
state that develops according to H a but is in 7i a} 

i^-Ha^j *„(a,<) = 0, (3.4) 

where the label a contains all the other information including the arrangement channel 
(even though including the arrangement channel in a is redundant for \I/ a (a,t), it is con- 
venient for other purposes). 

We then define *(+)(<*,<) state in 7i that develops according to H , 



(i^-Hj *(+> (<M) = 0, 



(3.5) 



and for which there exists an a-state such that 



lim U a (a,t),^(a,t)) = 1. (3.6) 



t^ — oo 



Therefore, the state \I/ a (a,t) is the "in" state \I/i n (a,t) in relation to the state ^^ + \a, t). 
Eq. (3.6) demands that the probability of finding the system in state ty a (a, t) in the remote 
past approach 1, and therefore, it is equivalent to 



* a (a,<), (3.7a) 



or 



j dt iG+ (t - t') *(+) (a, t') ^ (a, t) , (3.76) 



with the double arrow denoting the strong limit. 
Similarly, 

j dt' (-i) G~ (t - t') (a, t') *out (a, <) 



(3.8) 
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Exactly analogous to Eqs. (2.11), we can now write 

/+oo 
dt'G + (t-t')#> in (a,t') , (3.9a) 
-oo 

/+oo 
-oo 

We can now define the Moller matrices, and the S-matrix. The Moller matrices are 
defined by 

<f ± (a,t) = Q { a k h a (a,t), (3.10) 

with only those states \I/ a admitted which are in 7i a - On the orthogonal complement (i.e. 
the ionized eigenstates of H a ) fiW is defined to be zero, 

fi(±)p a = fi(±). 

Then, on the space 7i a} using Eq. (3.9a), we find 

(+) _ p , K (+) 

a1 +) = -i / dtg+(-t)H' a G- a (t)p a 

J — oo 


-i I dte lHt H' a e- lHat P a , 

oo 



and therefore, 



0l +) = lim e lHt e- lHat P a . (3.11) 

t^ — oo 

We can similarly find, on 7i a} that 

0l _) = lim e lHt e- lHat P a . (3.12) 

t^oo 

The range of fii + '' is the space of all full states that develop from arrangement channel 
a, and the range of fi„ ^ is the space of all full states that develop into arrangement channel 
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a. Let us call these ranges 1Z^ and izi \ and their respective orthogonal projections 
and Qa \ The Moller matrices, fiW, map 7i a onto TZ^, and from Eqs. (3.10) we find that 
on 1Z^ and TZa \ respectively, 

^ a (a, t) = ^ m (a, t) = (a, t) 

= ^out (a, t) = ftH V") (a, <) . (3.13) 

Therefore, because the ty a (a,t) span the space 7i a} we find that the Moller matrices, 
fil^, are partially isometric from the space 7i a , i.e. 



fii^fii^ = P a . (3.14) 



Similarly, the fil^ are partially isometric from the ranges ■]Z[ ±} of the fiW, i.e. 

n^W^ = Ql ±} , (3.15) 

which defines the The full states developing from or into any arrangement channel 

are orthogonal to each other as can be seen by direct evaluation of the inner products of 
asymptotic states. "If the two arrangement channels are different, then there must be at 
least one particle for which the "overlap" of the two states was negligible in the remote 
past because it belonged to a different fragment. Hence that inner product must vanish 
for all times" [16]. 

A major point of difference with our results from the Rearrangement Model is the 
statement, "note that the same argument shows that the inner product 

(tf 6 (/M),tf a (a,<)) (3-16) 

approaches zero as t — > ±oo (for a ^ b). But since H a ^ it is not independent 
of t and hence it does not generally vanish for finite times" (emphasis added) [16]. In 
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the Rearrangement Model, this is untrue: we show in section 7 that our states are all 
orthogonal to each other. 

From the Schrodinger equation, one can write 

#fii ±} =fii ±} # a , (3.17) 

which means that intertwines H and H a . This again is a major difference with the 
Rearrangement Model, because we show in section 7 that intertwines H and He, where 
He is the comparison Hamiltonian, which has the same spectrum as H; here, H a does not 
have the same spectrum as H . 

Our channel definitions could also include the single cluster arrangement channel, which 
is the channel of all the n-particle bound states of H . If we define A to be the orthogonal 
projection onto that subspace, then for all a we have 

Q^A = 0. (3.18) 

Now, every non-bound state must be decomposable into states that arise from, or go into, 
one of the other arrangements. Therefore, we assume 

A + ^<5i ±} = l, (3-19) 

a 

which is known as asymptotic completeness. 

Using Eqs. (3.12), (3.17), and (3.19), we may then define a unitary S-matrix, 

^ * * t (t) = Y,Sba*a(<X,t), (3.20) 

b 

where 

S ba =0(- }t 0i +) . (3.21) 
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The mathematical questions of convergence are the same here as for the single channel 
case, as in Eqs. (2.26). 



4. The Rearrangement Model 

To keep contact with earlier work, we shall use a combination of the notations of [7] 
and [17], as far as possible. We consider a quantum field theory with five distinct fields, 
B } C, -D, 6, and </>, and the corresponding particles (no anti-particles). 

The non-zero commutators are: 



i. 



v\<j>\v') =6{u'-u). 



(4.1) 



Note that 6 and <f> are labelled by continuum parameters, < u, v < oo, while B } C, 
and D } are treated as single modes ("infinitely heavy") [18]. We choose to use the energy 
as our variable, rather than momentum, because this makes the model much simpler, and 
more physically transparent. We want a total Hamiltonian for the system which allows 
the transitions 



B <-> C<j>, 



and 



d <-> ce. 



Therefore, we choose our Hamiltonian to be: 



H = H Q + V, 



(4.2) 
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where 



H = m B B^B + m D D^D + J du (u)6(u) + j ' dvv<j>\v)<j){v), (4.3) 



and 



V = J ' dto f {to)e{to)CD^ + Jduf* {to)e^{to)C^D 

+ J du g* {v )^{v)C^B. (4.4) 
This Hamiltonian has three constants of motion apart from itself. They are: 

d = B^B + tfc + D^D, (4.5a) 

<7 2 = B^B + J ' dv<f>i(v)<f>(v), (4.56) 

<7 3 = D^D + j du6i(u)6(u). (4.5c) 

Therefore, no transitions can occur between sectors labelled by different values of these 
quantum numbers. Let us start by enumerating the stable sectors. The first such sector 
is the vacuum and has C\ = C2 = C3 = 0. The next three are: C\ = 1,(72 = 0,(73 = 0; 
C\ = 0, C2 = 1, C3 = 0; and C\ = 0, C2 = 0, C3 = 1. These correspond to the states (7, <f>, 
and 6, respectively. Finally, there is the sector with C\ = 0, C2 = 1, C3 = 1; it corresponds 
to a state 6cf>. 

The three lowest non-trivial sectors are: 

(71 = 1,(72 = 1,(73 = 0, (4.6a) 

d = 1,<7 2 = 0,<7 3 = 1, (4.66) 

Ci = 1,C 2 = 1,C3 = 1. (4.6c) 

These correspond to B C<f), D C8, and B6 ^ C6cf> ^ D<f>, respectively. The last 
of these is the sector in which rearrangement collisions can take place, and as mentioned 
before, we shall call this the Rearrangement Sector. 
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Our strategy for solving the model in the Rearrangement Sector will be to first con- 
struct the solutions of the two lowest non-trivial sectors, (4.6a) and (4.6b) (which are 
exactly analogous to the Lee Model), and then use these solutions to express the Rear- 
rangement Sector equations. 



5. Solving the Model 

We start by constructing the solutions for the B <->■ Ccf> and D <->■ CO sectors. These 
are exactly the same as the Lee model, so the solutions are simple. We shall denote non- 
interacting ("bare") states by single bras and kets (( , )), and interacting states ("dressed" 
or "physical") by double bras and kets ((( , ))). 

The equations we need to solve for the continuum solutions are 

H\\))=\\\)) (5.1a) 



in the B <->■ Ccf> sector, and 



H\^))=^)) (5.16) 



in the D <->■ CO sector. 
We define 



a(z) = z-m D - I du * , (5.2) 

o z — u 

00 w 2 



f3 (z) = z — niB — / dv , (5-3) 

Jo z — v 

p X)B (v) = {C^v)\\)), (5.4) 

PlliD {u) = {C6{u)\ti), (5.5) 
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a XjB = {B\\)), (5.6) 
= (!%}}. (5.7) 

For shorthand, we will write a(X) for a(\ + ie) and a*(A) for a(\ — z'e), and similarly for 
/3(A). In terms of these, the solutions are: 

P\,b( v ) = b \ X ~ v ) + ~\ ~— ' ( 5 - 8a ) 

A — v + ze 

p„ >jD (a;) = * Qi - a;) + (5.86) 

= ^ (5 ' 8d) 

If a(z) develops zeros then we have additional discrete states. Similarly, if f3(z) develops 
zeros then again we have additional discrete states. For our purposes, we shall always 
assume that both a(z) and f3(z) have exactly one zero each, which are denoted by Mp 
and Mg, respectively. There is no loss of generality if we use this assumption because the 
extension to more than one zero is trivial. The equations for the discrete states are: 

H\M B )) =M B \M B )) (5.9a) 

in the B <->■ Ccf> sector, and 

H\M D )) =M D \M D )) (5.96) 

in the D <->■ C6 sector. We define 

P b(u) = (C<I>(u)\Mb)), (5.10) 

PD (uj) = (Ce(uj)\M D )), (5.11) 

y^=(B\M B )), (5.12) 

yfzE={D\M D )). (5.13) 



21 



In terms of these, the normalized solutions are: 



pb{v) = yZ B 
Z b 



M b - v 

r H 

M D - to' 



1+ dv 



\g{ v ) 



z D 



1 + / du 



l/MI 2 



(5.14a) 
(5.146) 

(5.14c) 
(5.14c/) 



(M D - uY 

where the last two are obtained by imposition of the orthonormality condition. Note that 
these solutions, Eqs. (5.8) and (5.14), form a complete orthonormal set. 

Now, we use these solutions to construct the solutions in the Rearrangement Sector. In 
this sector we will have four sorts of solutions. The first will correspond to the "physical" 
\C6{u)(f>(v))) sector, the second to the "physical" \D(f>(v))), the third to the "physical" 
\B8(uj))) } and the last to one or more dynamically generated bound states, which we shall 
denote by \M J s L )). Which solution is obtained will depend on whether we put the delta 
functions (which represent the plane wave parts of our solutions) in our solutions. If we 
put none, we get the discrete states. 

We wish to solve the eigenvalue equation 



H\E)) = E\E)). 



(5.15) 



We expand Eq. (5.15) in terms of 0t(u;)|A)), 8^ {u)\M B )) , ^t^)^)), and <ft{v)\M D )). By 
acting on these states with the Hamiltonian, and using Eqs. (5.1) and (5.9) we get 

(E - \ - u){{\\8(uYE)) = f H {{\\tfD\E)), (5.16a) 
(E — Mb — uj)((M b \8(ujYE)) = f* M ((M B \C^D\E)), (5.166) 
(E — [j, — v)(^\<t>(vYE)) = g* (v) ((HC<ti?|£}}, (5.16c) 
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(E-M D - v){{M D \^v)\E)) = g* (v) {{M D \tf B\E)) . 



(5.16c/) 



We need to evaluate the unknown matrix elements on the right hand side of Eqs. (5.16). 
We solve for these elements by inserting H in them, commuting it on one side, and letting 
it act on \E)) on the other. For example, we can solve for ((X\C^ D\E)) in the following 



manner: 



((X\C^DH\E)) = E((X\C^D\E)) 

((\\{HtfD + [ct£>,#] }|£}} = E((X\C^D\E)). 



(5.17) 



We now let H in the first term of Eq. (5.17) act on ((A|, and evaluate the commutator 
in the second term. We proceed similarly for the other three equations and, when all the 
dust settles, get 

(E- X-m D )((X\CD^\E)) 

du f (u) ((X\e(u)\E)) 



duf{u)[J dX'a xljB {{X'\e(u)\E)) + ^{{M B \e(u)\E))] 



+ J dvg(v)[J dv'a^,M\<l>(v)\E)) + yf^{{M D \<l>(v)\E))} 
(E-M B -m D )((M B \CD^\E)) 



(5.18a) 



dLof(Lo)((M B \6(u)\E)) 



duf(u)[j dX'o xlyB ((X'\6(u)\E)) + ^((M B \6(u)\E))} 



+ J dug(u)[j d^a^ >D ((^\<l>(u)\E)) + ^((M D \(l>(u)\E))} 
(E-fi-m B ){{fi\CB^\E)) 



(5.186) 



dvg(v){{^(v)\E)) 



duf(u) { J dX'* x , >B ((X'\6(u)\E)) + ^((M B \9(u)\E))} 
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+ j dug(u){j dfi'a^ D {{fi'\<l>(u)\E)) + ^{{M D \</>(u)\E))} 
(E-M D -m B )((M D \CB^\E)) 



(5.18c) 



dug(u)((M D \<f>(u)\E)) 



Z D 

+ 



duf(u)[j dX'o xlyB ((X'\6(u)\E)) + y^((M B \e(u)\E))} 



dug(u){j dfi'a^ >D {{fi'\<l>(u)\E)) + ^{{M D \<l>(u)\E))} • (5-18d) 

We first solve for the "physical" \C6{u)(f>(v))) states. We start by inverting Eqs. (5.16) 
and putting in the requisite delta functions. Note that we may put the product of two 
delta functions because this is an infinitely degenerate (double) continuum, which cannot 
be specified by just E; rather, we have to label the state with the variables E and n, with 
the n variable representing the division of energy between the 6 and <f> particles. We then 
substitute these into the first term of each of Eqs. (5.18), and solve for the unknown matrix 
elements. Having found them, we put them into Eqs. (5.16) to find our solutions. Defining 



b c \E,n,X,u) = ((X\e(u)\E,n)), 



b C F (E,n,M B ,uj) = ((M B \e(uj)\E,n)), 



d C {E,n,n,v) = (((j,\<f>(v)\E,n)), 



d c F (E,n,M D ,iy) = {{M D \</>(u)\E,n)), 



(5.19a) 
(5.196) 
(5.19c) 
(5.19c/) 



we get 



b C F {E,n,M B ,u) = 

d c (E, n, n,v) = 8{n - n) p E - n ,B{ v ) 



b (E } n } X } u) = S(E - X - n)p njD (u) 



(E - X - uj + ie) a (E - A) 



(E — M B — to + ie) a (E — M B ) 



K c {E,n), 



K c (E,n), (5.20a) 
(5.206) 



d$(E,n,M D ,u) 



(E - \i - v + ie) f3 (E - p) 



i /•: Mi, /' • /M ii i: m d ) 



K c (E,n). 



K c (E, n) , (5.20c) 
(5.20a 7 ) 
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where 



u (w \ 9(E-n)f(n) 1 

R c {E, n) = — — , (5.21) 

p [E — n) a (n) 7 [E) 

7 (E) = n/T f D „ , + fdfi \ f ^ ] }l * (5.22a) 

n ' (3(E-M D ) J ^ |a(/i)| 2 /3(£-/i) V ; 

- ZB A/A^ ! : ,' (5.226) 



a{E-M B ) J |/3 (A) | 2 a (E - A)' 
The last equality follows from Eqs. (A. 10). 

We now solve for the "physical" \Dcf>(v))) sector. We must again start by inverting 

Eqs. (5.16), but this time need to put just the one requisite delta function in Eq. (5.16d). 

We put these equations in Eqs. (5.18), and solve for the unknown matrix elements putting 

in another delta function in Eq. (5.18a) when inverting because now we must account for 

the zero of a(E — A) at Mj). We then put these results in Eqs. (5.16). Defining 

b D (E,X,u}) = ((X\e(u})\E)), (5.23a) 

b$(E,M B ,u) = {{M B \6(u)\E)), (5.236) 

d D (E,fi,u) = ((fi\(f>(u)\E)), (5.23c) 

d?(E,M D ,u) = ((M D \<f>(u)\E)), (5.23a 1 ) 

we get 

b°(E, A, u) = PD (c) S(E — A — M D ) - f ^ , ■ >. Z' B \ j <D (E) , (5.24a) 

(E — A — to + le) a (E — A) 



W E ' M °>^ = - { E-M B -, + >e )a{ E-M B ) KD<E) - < 5 ' 246) 



^E,M d , v) = pe _ MdMv) - [E _^ d ^ + k) ^J Md K d( E ) , (B.2M) 



where 



s/Z^ g(E-M D ) 
KD{E] = ^E)t>(E-M D Y (5 ' 251 
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and "y(E) is the same as that defined in Eqs. (5.22). 

In exactly the same way, we can find the solutions for the "physical" \B6(u))) sector. 
They are: 

b B (E, A, u) = - { {t °\., ^ B \ kb (E) , (5.26a) 
( ±j — A — to + ie) a ( ±j — A) 



AE^) = P B MHE-»-M B] - (E _^l + itji ^^. K B (E ] , (5.26.) 

where 

Kb(E) - 1 (E)a(E-M B y (5 - 2?) 
and ~f(E) is the same as that defined in Eq. (5.22). 

Finally, we wish to solve for any dynamically generated discrete states. In this case, 

we put no delta functions anywhere. When we follow the procedure of putting Eqs. (5.16) 

in Eqs. (5.18) and solving for the unknown matrix elements, we find that the only way 

to satisfy all the equations is if "y(E) has zeros. Denoting these zeros by M^, we find the 

discrete state solutions: 



im A ,M B ,.) = - {MA _ MB ^ + , e)a{] J A : Mg KA(M A] , (5.286) 



dfrU A , M D , „) = - {MA _ M ' B L v + ie)HM \_ MB) KA (M A ) , (5.28.0 

where Ka{Ma) 1S now an arbitrary normalization factor which is fixed, when demonstrating 



completeness, to be \J \e=M a ( see the discussion after Eq. (7.18)). For our purposes, 
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without loss of generality, we assume that there is only one zero of "y(E), denoted by M^, 
and thus only one dynamically generated discrete state. The extension to more than one 
discrete state is trivial. 

In each of Eqs. (5.20), (5.24), (5.26), and (5.28) the superscript refers to the sector that 
the solution is in, and the subscript F refers to solutions expanded in the the discrete state 
part of the Lee Model sectors. Furthermore, we have anticipated future developments by 
fixing the arbitrary constants accompanying the delta functions in Eqs. (5.20), (5.24), and 
(5.26). We do this by demanding that Eq. (6.3a) and Eq. (6.3b), or their equivalents for 
the other two sectors, give the same result, and that the solutions be orthonormal. 

6. Verification of the Solutions 

We now proceed to verify that Eqs. (5.20), (5.24), (5.26), and (5.28) are each solutions 
to our problem. To do this, we first transform our solutions into the bare state basis; 
i.e. in terms of the non-interacting states \C6{u)(f>(v)) , \B8(uj)) } and \D(f>(v)), using the 
completeness of the lower sector solutions. With the expansion coefficients in the "physical" 
\C6{u)(f>(v))) sector defined in the following manner (with the coefficients for the other 
sectors defined similarly) 



E,n)) = C c (E,n,Lo,v)\C6(Lo)^(v)) 




(6.1) 



where 




C C (E, n,u, v) 
D c (E,n,u) 
B c (E,n,to) 




(6.2) 



and 
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C C {E,n,u,v) = {C8(u)<f>(u)\E,n)) 
B c (E,n,u) = {B8(u)\E } n)) 
D c (E,n,v) = {D<t>(v)\E,n)). 



we have 



C C {E,n,u,v) = I d\px,B(v)b C (E,n,\,u) + pB(v)b%(E,n,MB,u), (6.3a) 



d^Pfi,D(^)d (E,n,p,v) + p D (u)d F (E,n,M D ,v), (6.36) 

c I c / — c 

D (E,n,u) = / dpa^^d (E ',n, p,u) + \J Zj)d F (E ',n,Mj) ,v), (6.3c) 

B c (E,n,u) = j d\<T\ tB b C (E,n, X,u) + \jz~B~b C F {E ,n,M B ,u), (6.3d) 

with similar equations for the other three sectors (for example, for the "physical" \Dcf>(v))) 
sector, we would replace C C (E } n,u,v) by C D (E,u,v), b c (E, n, A,cu) by b D (E, A,cu), etc.). 
A good check that we have solved our equations correctly is to verify that Eqs. (6.3a) and 
(6.3b) give the same result. This is indeed completely trivial if we use Eq. (A. 14). 
For the "physical" \C6{u)(f>(v))) sector in the bare basis, we get: 

C C (E, n,u>,v) = p niD {u)p E -n,B{y) ~ K C (E, n) f * ^ ^ M 



(E — to — v + ie) 
I^a,^! 2 , f j Wh,d\ 2 



dX — — — — + / dp 



(E - A - uj + ie)a (E - A) J (E — p — v + ie)/3 (E — p) 

+ a(E- Mb) (E — Mb - uj + ie) + (3 (E - M D )(E - M D -u + ie)^ 

D C (E, n, v) = f -^\p n , D (uj) - K C (E, n) g* (v) 
a ( n) 

d Wv,d\ 2 , i /g 4 ^ 

^ (E - p- v + ie)f3(E - p) (3{E - M D )(E - M D -v + ie) '' l ' ' 

B C (E, n, uj) = l {E ~ n \ pE-n,B(v) ~ Kc (E, n) /* H 
p [E — n) 

\ a \,B? | l (QA ) 

(E — X — uj + ie)a (E — A) a(E - M B ){E - M B - uj + ie) 1 ' C) 
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For the "physical" \Dcf>(v))) sector in the bare basis, we get 

n f * { '' 

C u (E,u,v) = p D {u)pe-M d ,b{v) ~ K D {E) 



(E — to — v + ie) 



dX — L — — — — + / dp 



(E — X — to + ie)a (E — X) J (E — p — v + ie)/3 (E — p) 

+ a(E- Mb) (E - M B - lo + ie) + /3(E-M D )(E - M D -u + ie)^ 
D D (E, v) = ^ PE -M D ,B (v) ~ K D (E) g* (v) 

dfi (E — p — v + ie)(3 (E — p) + (3 (E - M D ) (E - M D - v + ie) ^' ^ 6 ' 5 ^ 
B D (E, to) = p D (to) a E _M D ,B ~ K D (E) f* (lo) 

d ^ \ a \,B? , i (q k \ 

(E-X-UJ + ie)a (E-X) a(E- M B ) (E - M B - lo + ie) ' ' 1 ' C) 

For the "physical" \B6(lo))) sector in the bare basis, we get: 

C B (E,u, v) = PB {y) pe-m b ,d M - K B (E) f * {u) 9 * {v) 



(E — to — v + ie) 
dX—— ■ — ^^7— — — — + [dp ^ 



(E - X - lo + ie)a (E - X) J (E — p — v + ie)/3 (E — p) 

+ a(E — Mb) (E - M B -io + ie) + f3(E-M D )(E - M D -v + ie)}' 
D B (E, v) = p B (v) cr E -M B ,D ~ K B (E) g* (u) 

d/Jj (E — p — v + ie)(3 (E — p) + /3 (E - M D ) (E - M D - v + ie) } ' ^ 6 ' 6 ^ 

B B (E, lo) = ^p E -M B ,D M - K B (E) f* M 

dX Wx,b\ 2 Zb 1 , . 

(E — X — lo + ie)a (E — X) a (E - M B ) (E - M B - u + ie) J ' ' C) 

Finally, for the discrete states, we get: 



M^ — lo — v 
a X,B\ 2 , f j Wh,d\ 2 



dX ^ — -j— — + / dp 

+ 



(M A - X - to)a (M A - X) J r (M A - p-v)/3(M A - p) 
Z B 



a (M A - M B ) (M a Ma .o • /( i 
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+ /3 (M A - M D ) (M A -M D -u + ie) |' (6 ' 7a) 



D A (M A ,v) = -K A (M A )g*(v) 

{J dfi (M A - fi -v)f3(M A -fi) + P(M A - M D ) (M A -M D -u + ie)} ' ^ 
B A (M A ,u) = -K A (M A )f* (uj) 

f f dX W\,B\ 2 | ] (Q>j \ 

\] (M A — A — u)a (M A — A) a(M A - M B ){M A - M B - lo + ie) j' 1 ' C) 

We now verify that Eqs. (6.4), (6.5), (6.6), and (6.7) are each solutions of our model. To 
do this, we explicitly write down the analogues of Eqs. (5.16) in the bare basis, plug in 
each set of solutions in turn, and verify that the equations are satisfied. A straightforward 
analysis shows that the following equations must be satisfied in the bare basis (we have 
written them for the "physical" \C6{u)(f>(v))) sector, i.e. with the variable n — for the other 
sectors the variable n is, of course, missing): 

(E — u — v)C (E, n,u, v) = g* (u) B (E, n, u) + /* (u) D (E, n, v) , (6.8a) 
(E-m B -u)B(E,n,u)) = j dv g (u) C (E,n,u,v) , (6.86) 
(E-m D -v)D(E,n,v)= / du f (u) C (E, n, u, v) . (6.8c) 



Putting each of Eqs. (6.4), (6.5), (6.6), and (6.7) in turn into Eqs. (6.8), or their equivalents 
for the other sectors, and using Eq. (A. 14), we find that each of these sets of solutions 
satisfies the equations. Incidentally, a glance at Eqs. (6.8) immediately shows why we 
could not have solved the problem directly rather than the somewhat convoluted method 
we went through: the integral equations are not separable, and are quite intractable. 
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7. Orthonormality and Completeness 



We now proceed to verify orthonormality and completeness of the solutions Eqs. (6.4), 
(6.5), (6.6), and (6.7). We start by verifying orthonormality for the diagonal components 



beginning with the scalar product (\I/ C ' 1 (E' 



,n ,iu,z/] 



), ty c (E, n, u, i/)), which is given by 



dto dv ^ G 1 (E\ ri,u, u)V u (E, n, u, v) 



dcodvC u {E',n',u},v)C^{E, n, cu, v) 
+ j dwB c *(E',ri,u)B c (E,n,u) 
+ J duD c *(E',n',u)D c (E,n,u). 



(7.1) 



We now use Eqs. (6.3) to write this as 



dto dv ^/ c \e\ n',u, u)^ c (E, n,u, v) 



did dv 



dX' p\, B (v)b C \E\ n', X',u) + p%(v)b%' (E', n\ M B ,u) 



+ du 



+ 



dX px, B (v)b C (E,n,X,Lo) + p B (u)bp(E,n, M B ,to) 
J dX' a\, B b c \E\ ri, X',u) + yfZB~b c F \E', ri, M B ,u) 
j dXa XjB b C {E,n,X,to) + ^/z^b^{E,n,M B ,u}) 
dvD C *(E' } n' } v)D C (E } n } v). 



(7.2) 



We then do the integrals over A and A' to find 



dto dv V c \e', n',u, v)^ C (E, n,u, v) 



dXdwb c (E' } n\X } uj)b c (E } n } X } uj) 
+ J dub%\E\ri,M B ,u)b%(E,n,M B ,u) 
+ / duD c *(E',n',u)D c (E,n,u). 



(7.3) 
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Defining 

f*(n')g*(E'-n') 



Li (E',ri) 

, „ f (n) q (E — n) , 
L 2 (E,n)= J {, (7.4) 



a* (ri) f3* (E' - n'Y 

{n)g(E - n) 
a(n)/3(E - n)'' 



we find that the sum of the first two integrals is 

6 (E - E') S(n-n')-6(E'-n'-E + n) f * [ n ')J 

+ L\ (E\ n') L 2 (E, n) { ^ {EI) a ^ E , _ E + n) + 1 {E) a(E - E> + n>) 

L l (E'M)L 2 (E,n) \ -Z B f \*x, B \ 2 

>y*(E')j(E) \a*(E' - M B )a(E - M B ) J a* (E' - X) a (E - X) ^ ' ) 

while the third integral gives 

r(n')f(n) 



6(E' -ri -E + n) 



or 



(n ! ) a (n) 



+ 



Ll (E\ n>) L 2 (E, n) { ^ {EI) a ^ E , _ E + n) + ^ (£?) a (E — E' + n') 
L^E'.n^^iE.n) 



T(Ei)l(E) 



fi (E - M D ) a* (E' - E + M D ) 0* (E' - M D ) a (E - E' + M D ) 
+ / 8- (V - - V + „>) + / \ 



I 2 r \„ |2 



.(7.6) 



fi (E - fi) a* (E' - E - fi) 

Adding Eqs. (7.5) and (7.6) together, and doing the integrals by combining them into a 
single contour integral (which evaluates simply to its residues), we find that the only term 
left is S(E' — E)6(n l — n), which is just as required. 
We can similarly show that 



j du dvV D \E' ,u,v)V D (E,u,v) 

dX dto b D *{E', A, u})b D {E, A, to) 
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and 



+ J dubf(E',M B ,u)bf(E,M B ,u) 
+ J duD D *(E',u)D D (E,u) 

= 6{E'-E), (7.7) 

J dLodv^ B \E\Lo,v)^ B {E,Lo,v) 

= j d\du>b B *(E',\,u)b B (E,\,u) 
+ J dubf(E',M B ,u)bf(E,M B ,u) 
+ J dvD B *{E',v)D B {E,v) 

= S(E'-E). (7.8) 



Finally, 



* lt i.W !.„•./■< it I...'./'! 1. (7.9) 
Now we take up the off-diagonal elements. For 

J du dvV c \E' ,ri ,u,v)V D (E,u,v) 

= j d\du>b c *(E',ri,\,u)b D (E,\,u>) 
+ J dubp*(E',n' ,M B ,u)bp(E,M B ,u) 

+ j dvD c *{E',n',v)D D {E,v), (7.10) 

we find that the third integral exactly cancels the sum of the first two, giving us 0. We 
can similarly show that 

j du)dvV c \E',ri,u,v)V B (E,u,v) = 0, (7.11) 

dto dv^ c \E' ,ri ,u},v)^ A {M A ,u},v) = 0, (7.12) 
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J dLodv^ D \E\Lo,v)^ B {E,Lo,v) = 0, (7.13) 
J dLodu^ D \E',Lo,u)^ A (M A ,Lo,u) = 0, (7.14) 
J dto dv^ B \E' ,Lo,v)^ A {M A ,Lo,v) = 0. (7.15) 

Therefore, the set of solutions we found, Eqs. (6.4), (6.5), (6.6), and (6.7) are orthonor- 
mal. 

We now consider completeness. We wish to show that 

J dEdnV c (E,n,u),u)ty c \E,n,u)',u') + j dE ty D (E, to, v)^ d \e, J, v ) 

+ J dE^ B {E,u,v)^ B \E,J ,v') + ^ a {M a ,u},v)^! a \m a ,J ,v') 

' 6{ v -v' )6{to - to') \ 

8(u-u') . (7.16) 

6{lo-J)) 

Let us start with the diagonal elements. The (1, 1) element of the matrix is 

dE dnC c (E,n,u,v)C c \E,n,u',v') + j dEC D {E,u},v)C D *{E,J,v') 
+ j dEC B (E,u},u)C B *(E,u}\u') + C A (M A ,u},u)C A *(M A ,u}',u'). (7.17) 

These integrals are most easily done in the following manner. The first term can be 
rewritten, using Eqs. (6.3), as 

dEdn C C (E, n,u, u)C c *(E, n, J, v) 

d\' p\, B (y')b c \E, n, X',u') + p* B (y')b c F \E, n, M B ,J) 



dE dn I I dX p\ ! B( iy )b C '(E,n, + p B (v)b ( p(E,n, M B ,u) 



(7.18) 



One then rewrites subsequent terms in Eq. (7.17) in a similar fashion as Eq. (7.18). 
Since the integrals are exceedingly tedious, we describe how they are done, and leave it to 
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the interested reader to verify our results. The integrals over n are done with the help of 
Eq. (A. 15). Then, the integrals over E are done by converting them into contour integrals. 
When all the contour integrals are combined it is found that they add together into one 
large contour integral (plus the non-contributing circle at infinity), which evaluates simply 
to its residues. These residues exactly cancel the other pieces in the expression, leaving 
over one or more delta functions for the diagonal terms, and nothing for the off-diagonal 
ones. For convenience, the branch cuts and poles of the function ;^~y, where z is a complex 
integration variable in the contour integral, are shown in Figure 3. 

One finds that the (1,1) term is 8{u — u r )6(v — v r ). In doing this, one has to fix 
K A (M A ) = J dl }P \e=M a , which fixes the unknown normalization constant in Eqs. (5.28). 
One can similarly show that the (2,2) and the (3,3) terms are 6(v — v r ) and 8{u — u/), 
respectively. 

For the off-diagonal terms, one proceeds similarly and finds that they are all zero. Thus, 
our set of solutions, namely, Eqs. (6.4), (6.5), (6.6), and (6.7) is a complete orthonormal 
set of solutions of our model in this sector. 



8. The Moller Matrix and the Comparison Hamiltonian 

The matrix (with continuous eigenvalues) of the eigenfunctions, including any discrete 
solutions, gives us the generalized Moller Matrix by virtue of the results already demon- 
strated on orthonormality and completeness [7]. It is given by 

Sl(E,n,u,,u) = (^ C (E,n,u,v)^ D (E,u,v)^ B (E,u,v)^\M A ,u,v)) (8.1) 
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with components 



C c \E,n,u,v) C°i /•.".„•./.<! C"i /•.".„•./.<! C A {M A ,u,v) 
D c (E,n,u) D D (E,u) D B (E,u) D A (M A ,u) 
B c {E,n,J) B D {E,J) B B {E,J) B A {M A ,J) 



It has the properties of being unitary 



fifit = 
fit ft = 1, 



and of diagonalizing the full Hamiltonian, H , 



#fi = fiffc, 
fitfffi = 



(8.2) 



(8.3a) 
(8.36) 



(8.4a) 
(8.46) 



where He is called the comparison Hamiltonian. It can be calculated in the following 
manner. First, we use the eigenvalue equations to write 



HQ, (E } n, u, v) 



/•.Vi /•:. //.„•. r i. im> ij > /•:.„•. r i. /:t "i /•:.„•. /< .\/ 'i .\/ /< \ 



;D, 



(8.5) 

and then act on Eq. (8.5) with fit from the left, and make use of the orthonormality 
relations to get 



fitfffi 



/ES(E - E')S(n-n') 









\ 

ES(E-E') 

ES(E-E') 

M A J 



He- 



(8.6) 



To compare Hq with the free Hamiltonian, Hq } we rewrite Hq } putting E = n + r 
for the (1, 1) element, E = Mp + r for the (2, 2) element, and E = Mb + r for the (3, 3) 



(8.7) 
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element, and similarly for E' . Thus, He becomes 

/{n + T)6{T -T')8{n-ri) \ 

(M D +t)6(t -t') 

(M b + t)S(t -t') 

\ Ma / 
The free Hamiltonian, Ho, is 

' {to + v)6{lo -lo')6{v -v') \ 

(m D + u)8(u -v v ) . (8.8) 

(m B + lo)S{lo -J)) 

Comparing He and Hq } we see that we can identify He with Hq if we include both 



mass and wave-function renormalization terms in the interaction, and ignore the discrete 
Ma state in He- The mass renormalization means that we must add a quantity A to Ho, 
where A is 

/0 \ 

A= (M D -m D )6{v -v r ) . (8.9) 

\0 (M B - m B )6{io - J) J 

The structure of our solutions, Eqs. (6.4), (6.5), and (6.6) immediately tells us that we 
must have a wave function (and consequent coupling constant) renormalization. 
Thus, the fields B, C, -D, 6, and <J> have the wave function renormalizations 



B -> JffB = -^=B, (8.10a) 

D -»■ V^D = —^=D, (8.106) 

C -> C, (8.10c) 

-> 0, (8.10c/) 

<^-><^. (8.10e) 

Because there are no proper vertex corrections, the coupling constant renormalizations 
reflect the wave function renormalizations [7] 

/( W )^V^/( W ), (8.11) 

g(v)^yfZ^g(v). (8.12) 
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Furthermore, as there are no divergences in this problem, the coupling constant and 
wave function renormalizations are inessential, and the mass renormalization making Hq + 
A identifiable with Hq is essential only in this sector. These renormalizations are sufficient 
for higher sectors as well. The only change in the higher sectors is due to the mass 
renormalizations which alter the continuum thresholds from mp and mjj to Mp and Mg, 
respectively, but leave everything else unaffected. 

Notice that while Hq and Hq have the same structure (as long as a and f3 both 
have zeros, and 7 does not), they have different spectra. Only the double continuum 
0<n<S<oois coextensive; the Dcf> and B6 continua are renormalized downwards from 
mj) to Mj) and from mjj to Mg, respectively. Notice also that, contrary to conventional 
wisdom [1,3,12], the Moller matrix intertwines the full Hamiltonian, H , with Hq } not with 
Hq. However, Hq and H do have the same spectrum. 

In addition, if we take the unitary transformation of He in reverse, we can convert the 
comparison Hamiltonian to the full Hamiltonian 



and just as in the Cascade Model of [7], we find that the notion of an interaction is basis 
dependent. 




(8.13) 
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9. The S Matrix 

We have obtained one set of solutions to our problem, namely, Eqs. (6.4), (6.5), (6.6), 
and (6.7). We can, of course, obtain another set in which the singular operators of the form 
(E — to — v + ze) _1 (which give the "in" states) in Eqs. (6.4), (6.5), and (6.6) are changed 
to (E — to — v — ze) _1 (which give the "out" states), while Eqs. (6.7) remain unchanged. 
Let us denote these solutions, and quantities associated with them, with a prime. This 
new set also furnishes a Moller matrix, 



which satisfies the same properties as the original Moller matrix, that of unitarity 



0' = (^! C \E,n,u,v)^ D \E,u,v)^ B \E,u,v)^ A (M A ,u,v)) , 



(9.1) 




= 1, 



(9.2a) 




= 1, 



(9.26) 



and of diagonalizing H 



HQ' = Q'H C , 



(9.3a) 




He- 



(9.36) 



The set of states, Eqs. (6.4), (6.5), (6.6), and (6.7) are such that 

/ ' 6{n-Lo)6{E -to -v) 
lim e lHct e- lHt y C (E,n,LO,v) = 




(9.4a) 



(9.46) 



(9.4c) 



lim e ltict e 



*H c t „-iHt,-r,A 



\yfZ^6{E - M B -oj) 



(9Ad) 
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of which the first three are the plane wave ideal eigenstates of the comparison Hamiltonian. 
However, notice that there is the need for a wave function renormalization in ty D (E,u, v) 
and ^ B (E 1 to 1 v\ and that the threshold is renormalized in these two cases (i.e. m B — >■ M B 
and m B — > M B ). Clearly, these states are the "in" states in our problem. This is again 
analogous to the Cascade Model of [7]. 

For t — > +00 for these "in" states we have 



lim 



( 



6{E 



to 



y U U > a(n)P(E-n) + ^ ' ' " ' ' ' " ' 



2gi_ /(n) g (£-n) _ Md 



~/(E) a(n)P(E-n) a(w)P(y) 

9*(y) 



\ 



~/(E) a(n)P(E-n 
2-ki f(n)g(E-n 
-/(E) a(n)P(E-n 



^HMM^lz B S(E-M B 



a- 

a(to) 



lim 



e lHct e- lHt ^ D (E,uj,v) 



to 



( 2m8{E 
s/Z^S (E - M D 



,.\ VZ^9(E-M D ) f*(w)g*(E-w) 

> ~/(E) P(E-M D ) a(uo)P(E-uo) 



\ 



') 



P*{y) 



+ 



2-ki 



Z D - 



\ 



r^8{E-M B -u)$jkZ B 



l(E) " a(io)a(io) 



lim 



e lHci e- lHi ^ B {E^,v) 



( 2ttiS(E 



u — V 



VZ^f(E-M B ) f*( u )g*(E-w) 
-/(E) a(E-M B ) a(to)P(E-Lo) 

2wi | g(t/)| 2 



\ 



fZ D 8(E-M D -v)^Z D 



y(E)""P( V )P( V ) 



^y/Z~B~8(E -M b -uj) 



a* (to) , 2-ki 



+ 



-Z B 



lim 



a(io) 1 l(E) *~ ' 13 a(io)a(io) 



\ 



(9.5a) 



(9.56) 



(9.5c) 
(9.5c/) 



(The limits in Eqs. (9.4) and Eqs. (9.5) are understood for multiplication by smooth 
functions of u or v or both, as the case may be). 

The "out" states behave in an analogous but opposite fashion to the "in" states. They 
behave simply for t — > +00, but have a complicated structure as i -> —00. Furthermore, 
the "in" states at t — > —00 and the "out" states at t — > +00 are identical. Therefore, we 
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can define an S-matrix, and can compute it in one of several ways. For example, we can 



compute it using 



^scattered = Hm (t) - * (-<)) 



(9.6) 



or we can take the scalar product of the "in" and "out" states 



(9.7) 



Both methods, of course, give the same answer. 

The method of the inner products is cleaner and more aesthetically satisfying so we 
shall follow it for the calculation. The results are easily checked by doing the calculation 
by the other methods. 

Schematically, the S-matrix looks like (the "+" subscript means an "in" state and the 
"— " subscript means an "out" state) 

/-((ce<f>\ce<f>)) + -((ce<p\D<p)) + _{{ce<f>\B6)) + -((ce<p\M A ))\ 



s 



((D<P\Ce<P)) + _{{D<f>\D<f>)) + -((D<f>\B6)) + -((D<P\Ma)) 

((Be\ce<f>)) + -((Be\D<f>)) + -{{Be\Be)) + -{{bo\m a )) 

V ((M A \Ce<P)) + ((M A \D<P)) + {{M A \B6)) + ((M A \M A )) J 



Let us start with the (1,1) component of S. We wish to calculate 



(9.8) 



c((E' , n', out\E, n, in)) c 



n, u, v) 



+ J dioB c '*(E',n',io)B c (E,n,i l j) 
+ / duD c '*(E',n',u)D c (E,n,u). 



(9.9) 



We rewrite Eq. (9.9) in terms of the lower sector physical states using Eqs. (6.3) to get 



did dv 



dX 1 p'\, B {u)b c, \E' , n',X',u) + p'sWbC'/iE', n', M b ,lo) 



41 



+ / dto 



dX p\, B (v)b c (E, n, A, to) + p B (u)bp(E, n, M B ,to) 



dX' a' xljB b c '\E', ri, X',u) + ^Z^b c ' F \E', ri, M B ,to) 
j dX (J XjB b C (E, n, A, to) + y^b^E, n, M B ,u) 
+ J duD c '*(E',ri,u)D c (E,n,u). 

Doing the integrals over A in Eq. (9.10) we get 



(9.10) 



dXdX' ^-^-S (X- X') f dwb c '*(E',ri,X',u))b c (E,n,X,u)) 
P (A) J 

+ J dcob C 'p \E',ri,M B ,u)b%(E,n,M B ,u) 
+ / dv D C '*(E', ri , u)D c (E, n, v). 



(9.11) 



The sum of the first and second integrals gives 



fi(F-F'Mf>( - >) (E - n) a* (n) 2vrz g (E - n) f (n) g* (E - ri) /* (ri) 
K J ' 1 ' f3(E-n)a(n) 7 (E) f3 (E - n) a (ri) f3 (E - ri) a (ri) 

^(E-n)t(n>)f(n) 



+ 



g*(E'-ri)f* (n')g(E-n)f(n) 
(3 (E 1 - ri) a (ri) (3(E-n)a (n) 



+ 



+ 



1 



7 (E) a (E - E' + ri) 7 (E 1 ) a (E 1 - E + ri) 
g*(E> -n')r{ri)g{E-ri)f{ri) 



7 (E 1 ) 7 (E) (3 (E 1 - ri) a (ri) (3(E-n)a (n) 

1 



dX \a\ - |2 



Z B 



,5| 



a (E 1 - X) a (E - X) a (E 1 — M B ) a(E — M B ) 



(9.12) 



and the third integral gives 



6 (E' - ri - E + n 



P*(E-n)r(ri)f(n) 

f3 (E — n) a (ri) a (n) 
g*(E'-ri)f* (n')g(E-n)f(n) 
(3 (E' - ri) a (ri) (3(E-n)a (n) 



+ 



1 (E)a(E - E' + ri) 7 (E 1 ) a (E 1 — E + n) 
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1 g*(E'-n')f*(n')g(E-n)f(n) 
7 (E 1 ) 7 (E) (3 (E' - ri) a (ri) (3(E-n)a (n) 



+ 



E + fi) a(E'-E + (j.) 
1 f , 2 1 / 1 1 

+ 2j 11 (3(E'-ri) \a*(E-E' + ri) + a(E-E' + ri) 

Zp ( 1 1 



2(3(E-Mp) \a*(E'-E + Mp) a(E' - E + Mp) 
2f3(E'-Mp) \a*(E-E' + M D ) a (E - E' + Mp) J f ' y ' ' 

Adding Eqs. (9.12) and (9.13), and converting the sum of the integrals to contour 
integrals (which evaluate to their residues and cancel the other terms with them inside the 
curly brackets), we are left with 

c((E', n', in\E, n, out)) c = 6 (E - E') [s (n - ri) ^ ~_ 

2vrz g(E-n)f(n)g*(E-ri)t(ri) \ 
- f (E)/3(E-n)a(n) f3 (E - n') a (n 1 ) j' l ' ' 

In a similar fashion, we can do all the other S-matrix elements. They are 

f3* (E — Mp) 



D {{E',out\E,in)) D = 8{E-E')^ 



fi(E-Mp) 



B ((E',out\E,in)) B = 8(E-E'){ 



2iriZp |ff (E — Mp) | 2 
^ 1 (E) P(E-M D )P(E-M D ) h { ■ } 
a* (E - M B ) 



a(E- M B ) 



A 



+ 2 ^B \f(E-M B )? 1 

j(E) a(E - M B )a(E - M B ) )' 

((M A \M A )) A = 1, (9.17) 



I IT?' +\ TP ■ \\ 9 'A ( tp' ^ y/ZBVZ^ f* (E - M B )g(E - M D ) 
B ((E, aut\E, l n)) D = 2mS (E - E) a{E _ M]}) p {E _ M]) y (^) 

//tp' +\ tp ■ \\ 9 ^ y/Z^y/Z^f{E-M B )g*{E-M D ) 

D ((E, aut\E, m)) 5 = 2mS (E - E) a{E _ M]}) p {E _ M]) y (^) 

1 1 TP' +\tp ■ w o jp\ 9 (rc) f(E-n) g* (E - Mp) 

D ((E, aut\E, n, m)) c = 2mS (E - E) ^ ^ ^ _ r ^ _ ^ , (9.20) 
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c((E', ri, out\E, in)) D = 2tti6 (E 1 - E) 
B {{E', out\E, n, m)) c = 2kiS (E 1 - E) 
c{{E', ri, out\E, m)) B = 2ttiS (E' - E) 
A {{M A \E,n,m)) c = 0, 
c ((E',ri,out\M A )) A = 0, 
A ((M A \E,in)) D = 0, 
D {{E',out\M A )) A = 0, 
A ((M A \E,m)) B = 0, 
B ((E',out\M A )) A = 0. 



y/Z^g*{ri)r{E-ri)g{E-M D ) 
j(E) a(ri)(3(E -ri) (3(E-M D )'' 
y/Z^g(ri)f(E-ri) f* (E — Mb) 
7 (E) a (n) (3( E-n ) a( E- M B ) ' 

JZb~ g* (ri) /* (E - ri) f (E - M B ) 
j(E) a(ri)(3(E -ri) a(E-M B y 



(9.21) 
(9.22) 

(9.23) 
(9.24) 
(9.25) 
(9.26) 
(9.27) 
(9.28) 
(9.29) 



10. Unitarity of the S Matrix 

We can almost trivially show that the S matrix that we have obtained is unitary. In 
equations, we wish to show that 

SS^ = 1. (10.1) 
Let us calculate the (1,1) term in It is 



cot 

(U) 



6 (E — E') I dn" 



n — n 



+ 



8 (ri 



//A /3* (E — n) a* (n) 
P(E- n)a(n) 
2ttz r(n")f(n)g*(E-n")g(E-n) 
7 (E) a (ri 1 ) a (n) (3 (E - n") (3(E-n) 
f3(E -n')a(ri) 



n 



P* (E - ri) a* (ri) 
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(2mY 



;8(E-E 



2ttz f(n")f*(n')g(E-n")g*(E-n') 
7* (£) a* (n") a* (ri) (3* (E - n") (3* (E - ri) 

„ f(n)g(E-n) /* (ri) g* (E — ri) 
a (n) /3(E-n) a* (ri) /3* (E - ri) 



\g(E-M D )\ 2 \f(E-M B )\ 2 
-— + Zjb 



\/3(E-M D ) 



\a(E-M B )\ 2 



(10.2) 



Doing the integral in Eq. (10.2) with the help of Eq. (A. 15) we find that the result is 
S(E — E')S(n — n'), exactly as desired. The rest of the terms are done in the same way. 
We find 



SS} 2>2) =S(E- E>) 



9 — 1 



(10.3) 
(10.4) 
(10.5) 



with all other terms in SS^ being zero, as required. Thus SS^ = 1. In the same way, we 
also show that S^S= 1, and therefore, our S-matrix is unitary. 



can 



11. Eigenphases of the S Matrix 

The interesting case for the S Matrix is when E > so that all channels are open. The 
S Matrix must satisfy 

SC = rC, (11.1) 

where |r| 2 = 1, for some C. This is equivalent to the following relations (where we ignore 
the discrete A channel, as it is decoupled from everything else, and suppress S(E — E')) 

P* (E -n)a(n) . 
T ~ P(E-n)a(n) 
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2m f(n)g(E - n) f f , f*(n')g*(E - n') , ^(E-Mp) , 
>y(E)a(n)0(E-n)\J U a(n')0(E - n') n V D 0(E — Mp) 



+ VZB a(E - M B ) QB (1L2a) 



P*(E-M D ) C 
T 8(E — Mp) 
2m 



l(E) 



— g(E-M D ) \ f , f*(n')g*(E-n') , g*(E-M D ) t 
Zd P(E-M d )\] dU a(n>ME-n>) ^ + ^ 0(E - M D ) C 



D 



^- f\E-M B ) , , 



T a(E-M D ) B 
2m 



l(E) 



— f(E-M B ) \ f , r{n')g\E-n') r ^- g\E - M D ) , 
a(E — M B ) \ J ^ a(n>)8(E-n>) ^' + 0(E - M D ) ^ D 



/-^— f*(E — Mb) f , 
+ VZB a(E-M B ) (By (1L2C) 



We now define the unimodular quantities 



, N 8 (E — n) a* (n) 

Tin) = —-7- ( // , (11.3a) 

V ; 0{E-n)a{n) ' 

8* {E — M D ) 

TD = HE-M D y 

_ a*(E-M B ) 

= a(E — Mb) (1L3C) 
These are the basic equations. We can solve them for continuum values or for dis- 
crete values of the eigenphase shifts. Let us start with the continuum values. We invert 
Eq. (11.2a) and put a delta function on the right hand side along with the appropriate 
normalization. We then multiply both sides of the equation by ^ J'^^E-n) anc ^ integrate 
over n to get 



\f(l)\ 2 \g(E-l)\ 2 r(Q } f dn , rW(E-n') c 



2m 

TfE) J \a(l)\ 2 \8(E -1)\ 2 t -r(/) + zeJ J a(n')0(E-ri) 
r(T)g*(T) + 2m f \f(l)\ 2 \g(E-l)\ 2 r (I) 



<x(t)0(t) 7 (E)J \a(l)\ 2 \0(E-l)\ 2 T-r(l) + ie 



g\E-M D ) c ^ r(E-M B ) c , 
Z ° 0(E - M D ) Q ° + VZB a(E - M B ) QB (1L4) 
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Defining 



2ttz \g(E-M D )\ 2 a D 
-Z D - 



7 (E) u \P{E-M d )\*g-(j d 



2vrz l/^-M^)! 2 o B 



7(E) |a(E - M 5 ) | 2 a - a B 
we invert Eq. (11.2b) and Eq. (11.2c) to solve for the term in the curly braces on the right 
hand side of Eq. (11.4), namely, V^jj§Ej^^D + VZb ^e-Mb) ^-' and find 



Z ° 0(E — Mb) + VZB a(E-M B )^ B 
I_ A f , , nn')g\E-n') r 
S 7 7 a(n')/3(S - n') 



(11.6) 



We now define 



1 j " 7 (^)S7 d >(/)| 2 |/3(^-0l 2 ^-^(0 + ^' 1 j 



and use this to combine Eq. (11.4) and Eq. (11.6), and find 

J dU a(n')0(E-n') ^' + VZ ° f3(E — Mp) + V ^ a(E - M B ) 
1 r{T)g*(T) 



B 



X(t)?: a (t) P (E - t)' 
Therefore, our continuum solutions are 



(11.8) 



C n = Vr'S(T-r(n)) 

2vrz f(n)g(E-n) 1 f* (t) g* (E - t) 

1 {E) a(n)0(E - n) (r - r(n) + ie) X(r)S a (r) (E — t) ' 1 J 
_ 2vrz ^- g (E-M D ) 1 r(r)gHE-r) 

D 7 (E) V i? /3(E-M i3 )(r-r i3 + z e )X(r)E a(r)/3(E-r) ' 1 ' J 

C _ 27rz f^fi E - M B) 1 /*(r)flf*(E-r) 

~~ t(E) 5 a(E-M 5 )(r-r 5 + ze)x(r)E a (r)/3(E-r) ' (H ' 9c) 

To investigate the spectrum of r, we use the method of [7]. We define the following 
quantities, taking advantage of their being unimodular: 



e 



2i6, W = r(n), (11.10a) 



47 



e 2l ° D = t d , (11.106) 

e 2l0B = t b , (11.10c) 

e 2lh = t. (ll.lOd) 

We then put these definitions in Eqs. (11.9), and see that r(n) = e 2l(> ( n ) ranges contin- 
uously along a unit circle in the complex plane from = 0(0) to = 0(E). 
In addition, these solutions are continuum normalized 

dn C n (t - ie) C n (t + ie) 

+ £ D ( T ' - « e ) £ D ( T + « e ) + ^B ( T ' - « e ) ^B ( T + « e ) 

= 6{t'-t), (11.11) 

and will be complete if there are no discrete zeros of X(r). If there are, they will have to be 
included in the completeness identity. We now find the number of discrete zeros of X(r), 
that is, the number of discrete eigenphase shifts of our S-Matrix. 
We define 

1 + ix 1 + ix (n) 

T = i T ( n ) = i ryr, 

1 — ix l — ix ( n) 

_ 1 + ix p _ 1 + ixfi 

T D = Z : , T B = ~ : , 

1 — ix j) 1 — ix b 

put these in Eq. (11.7), and take real and imaginary parts to get 

1 f dJ \f(l)\ 2 \g(E-l)\ 2 x(l + x(l)) 



2 J \a(l)\ 2 \f3(E-l)\ 2 x-x(l) + 



ie 



Im 7 (E)\ l x(l + xp) ^ \g (E - Mp) | 2 



2vrz J 2 x-x D \P{E-M D )\< 



l,(l +XBlz \f(E-Ms,f 



2 x-x B \a(E - M B )\ 2 ' 



l - fdl \fV)\ 2 \9(E-l)\ 2 



2 J \a(l)\ 2 \(i(E-l)\ 2 

-Ke( l{E) \ l 7 \9(E-Mp)\ 2 1 \f(E-M B )\ 2 
~ Re ) ~ 2 Zd \HE-M d )\ 2 - 2 ZB \a(E-M B )\ 2 - 
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We observe that Eq. (11.12b) is an identity, by means of Eq. (A. 15). To find the number 
of zeros of X(r), we multiply both sides of Eq. (11.12a) by (x — xj))(x — 15). We then 
find that the highest power of x appearing in Eq. (11.12a) is x 3 , barring any higher powers 
contributed by the integral. Therefore, there are at least three discrete zeros of X(r), and 
thus, at least three discrete solutions which will have to be included in the completeness 
identity, Eq. (11.11). 

12. A general formalism for scattering theory 

In this section, we describe an approach due to Sudarshan and collaborators [19-25], 
which takes a very different view of scattering problems, and is quite different in spirit. It 
is essentially immune to many of the problems that occur in the conventional approaches 
in the literature. The idea is that one always works with the complete set of eigenstates 
of the full Hamiltonian, H , properly labelled. This set, by definition, is both orthonormal 
and complete. The matrix made up of these eigenstates is the Moller matrix. This Moller 
matrix, again by definition, will diagonalize the full Hamiltonian giving the comparison 
Hamiltonian. 

In other words, we have a Hamiltonian for the system, H . We wish a physical inter- 
pretation of this object as a scattering system. One starts by considering the complete set 
of states for H , which we denote as ij) a (E) } where E is the energy of the eigenstate, and a 
contains everything else necessary to uniquely specify the state such as spin, channel, etc. 
Then, 

H^ a (E) = E a ^ a (E). (12.1) 
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Form the generalized Moller matrix, W, by defining 



W E , a = V>« (E) . 



Therefore, 



HW = WH C . 



(12.2) 



where the implied integration is of the Stielje's type-z.e. we sum over any discrete in- 
dices, and integrate over any continuous ones. Because we have assumed that the set of 
eigenstates of H is complete, this Moller matrix has the property that 



It is thus isometric and unitary, as long as we ensure that the spectra of H and Hq are the 
same, and the spectrum multiplicity is properly preserved. The only caveat here is that 
not all formally hermitian Hamiltonians have a complete set of eigenstates. However, all 
"reasonable" Hamiltonians will have a such complete set. 

Now, one normally wants an interpretation of a scattering system in terms of asymp- 
totic states. The point here is that, to get such an interpretation, we should use He, not 
Hq. We must set up a correspondence between the set of eigenstates of He with H , because 
unlike Hq } we are guaranteed that H and He are isospectral. Any asymptotic conditions 
(such as the strong convergence properties of Eqs. (2.26)) should be expressed using Hq } 
not Hq. Thus, it is Hq } not Hq } which is the proper starting point for any perturbative 
scheme. Furthermore, as indicated by our model, we should endeavour to construct a per- 
turbative scheme to calculate the full states, not the asymptotic ones, because we cannot 
say, with any confidence, what the appropriate conditions are on the asymptotic states 
(see the next section for details). In addition, note another advantage of this formalism. 



fifit = i 5 oto = i. 



(12.3) 



50 



Any, and all, shifts in the thresholds and spectrum of H (such as a mass renormalization) 
are automatically taken care of by this procedure. 

Finally, it is worth remarking that, in general, W and Hq are not going to be analytic 
in the coupling constants. Therefore, the procedure that one occasionally sees in the 
literature of splitting H into He + V is not very useful, and not very constructive. Both 
He and V' will be complicated functions of the coupling constants, and will have all sorts 
of renormalization factors appearing. 

If we are interested in working perturbatively, then we must set up the system carefully. 
We give here an analysis of Sudarshan [26]. Consider a quantum system defined in a Hilbert 
space 7i with the Hamiltonian split in the usual way 



in which we already know the ideal eigenstates for the continuum, and the proper eigen- 
vectors for the discrete states. The ideal states are, of course, not normalizable and we 
must take proper linear combinations of them to get states that are square integrable, and 
in 7i. Then, we set up a correspondence between eigenstates of H and eigenstates of Hq } 
in such a manner that 



H = H + V } 



(12.4) 



(12.5a) 



#0^0 A = A?/>oa- 



(12.56) 



Therefore, 



(1-Go(A)T/)Va = V>oa, 



(12.6a) 



G (\)(\-H ) = l, 



(12.66) 
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where Go is the free Green's function. From this we can write 

V'A = (l-Go(A)y)- 1 VoA, (12.7) 

thus defining for us a possible Moller matrix, fi' given by 

O' = (1 - GoVy 1 . (12.8) 

Now this fi' is a possible Moller matrix in the sense that it intertwines H and H$: 

HQ 1 = 0,'Ho. (12.9) 

Unfortunately, it is not very useful because it is not necessarily unitary, or even isometric. 
One must renormalize it correctly so as to get a unitary operator. Furthermore, H and 
Hq are not isospectral. Consider the full Green's function, £7(A): 

Q (A) = x _l + %e = (1 - G (A) V)- 1 Go (A) . (12.10) 

While, at first glance, it would seem that Go and Q have the same singularities, this is 
not necessarily true. Firstly, Q can have additional singularities from the first factor (1 — 
Go(A)y)- 1 inEq. (12.10). These can come from bound states produced by the interaction, 
and more importantly, from continuum states in which one or more of the particles is 
composite so that its mass gets shifted. Secondly, Q can have some of its singularities 
cancelled when this factor vanishes. Therefore, Q and Go are not necessarily isospectral, in 
general. In other words, the statement that "perturbations vanish at infinity" is not valid 
generally. Rather, this naive asymptotic condition is not generally fulfilled. This shows us 
why ft 1 failed to be unitary: the new spectra produced by (1 — Go(A)V) -1 do not appear 
with a canonical weight. As advertised, we correct this problem by defining a renormalized 
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O given by 



O = (1 - G V)~ L D 



-l 



(12.11a) 



where 



D 2 = (1-VG+) L (1-G V) 



-l 



(12.116) 



However, since this new O is unitary, it connects H with an isospectral and diagonal 
Hamiltonian. We have already seen that this associated diagonal Hamiltonian cannot be 
Hq. Rather, it is a different object, which we call the comparison Hamiltonian, Hq- For 
further details, and concrete examples of this formalism applied to several models, such as 
the Lee Model, the separable potential model, and the Cascade model, see [26]. In these 
models, one can explicitly see these various effects such as shifts in the continuous spectra, 
the deletion of spectra from Hq to get the spectra of H , and the augmentation of spectra 
in Hq to get the spectra of H . 

13. Putting the "generic" formalisms to the test 

There are many different approaches to Quantum Scattering in the literature. The most 
familiar of these is potential scattering. Others include the LSZ formalism, the "almost 
local" formalism, and the Lax-Phillips formalism. Lax-Phillips [27] theory is outside the 
scope of this work. 

The well known LSZ formalism [5], extended by Mohan [28], postulates the convergence 
of the matrix elements of interacting fields to the matrix elements of free fields. However, 
the formalism does not apply in many cases. For example, as noted by LSZ themselves, it 



53 

is inapplicable to problems in which stable bound states exist. Trouble occurs when this 
point is forgotten, and the formalism is extended into areas where it is inapplicable. The 
"almost local" formalism due to Haag [1], Ruelle [9], Ekstein [10], Jauch [12], Araki [29], 
and others tries to be general enough to consider complicated problems [1]. Its basic 
idea is that it is possible to construct asymptotic in-going and out-going states as strong 
limits in Hilbert space, if a certain "space like asymptotic condition" is verified by the 
vacuum expectation values of products of field operators [9]: the so called "almost local" 
operators [1]. 

We shall restrict our attention to the conventional, and quite "generic" formalism, 
as reviewed earlier in sections 2 and 3; as mentioned before, the LSZ formalism is not 
applicable to situations where stable bound states are present, such as our model. We will 
compare these results to the results obtained from the Rearrangement Model. 

Conventional formalisms for quantum scattering theory have the following protocol for 
generic scattering systems: 

1. They do not use the comparison Hamiltonian 

2. The asymptotic states are orthonormal [1,9,13] 

3. The completeness of the asymptotic states is postulated [1,13] 

4. For the case of potential scattering only, the Moller Matrix is isometric but not 
necessarily unitary [1,13] 

5. The eigenstates of the exact Hamiltonian are never considered. 

We shall take up these points one by one, and put them to the test by comparing them to 
the results explicitly obtained from our model. 



1. It is essential when taking the limits, lim.t^± 00 e lIdi e - lId oiq> ^ where ^> is either 
a wave function or a field operator, that the continuous spectra of H and Hq 
coincide. If they did not there would be wild oscillations while taking the 
limit, and the limit would not exist. It is for this purpose that Hq is mass- 
renormalized to Hq. However, in general, this is still not enough. It is perfectly 
possible, if there are bound states or unstable particles in the spectrum of H , 
that no amount of tinkering with Hq will make its spectrum coincide with 
H. This can be seen by inspection of Eqs. (8.7) and (8.8). No amount of 
renormalization of Hq can give us the discrete state present in He, but 
this may be ignored because Ma is a discrete point eigenvalue. On the other 
hand, we do have the possibility of a continuous spectrum in H corresponding 
to the scattering states involving physical B or D particles. 

However, unlike H and Hq } H and He are guaranteed to be isospectral 
because He is obtained by diagonalizing H. Therefore, it is He, and not 
Hq } that is the proper starting point for any scattering scheme, perturbative or 
otherwise. The method for obtaining the correct spectrum of H by perturbation 
theory is discussed in the work of Sudarshan, Chiu, and Bhamathi [30]. In 
simple cases such as when stable bound states are not present, or field theory 
with no bound states or unstable particles, He can be identified with the 
renormalized Ho, as noted in section 8. 

In fact, even in cases where (formally) no splitting is made, i.e. no explicit 
mention or use is made of an Hq, there is still the implicit use of Hq because, 



commonly, asymptotic particles are defined as solutions of free particle equa- 
tions like the Klein-Gordon equation. 

2. Both formalisms assert the orthonormality of the asymptotic states, and the 
result is supposed to be generic. In Eqs. (9.4), we have obtained the asymp- 
totic states of the Rearrangement Model according to both formalisms. Yet, 
as we can see from a glance at the Haag-Ruelle asymptotic wave functions, 
Eqs. (9.4), the asymptotic states computed according to their rules do not 
form an orthonormal set. This point should not cause confusion. Our full 
states, namely, Eqs. (6.4), (6.5), (6.6), and (6.7) are, indeed, all orthonormal to 
each other, as was shown in section VII. As a result, we have orthonormal sets 
of "in" and "out" states. However, when we calculate the Haag-Ruelle type 
asymptotic states according to either of the formalisms, we find that they are 
not orthonormal. This lack of orthonormality stems from a factor of the wave 
function renormalization constant that appears in each of the asymptotic wave 
functions. This factor is essential: if it were not present, the interacting states 
would not be orthonormal. 

3. As mentioned earlier, Ruelle extends Haag's work by postulating the complete- 
ness of the "in" and "out" states [9]. This is also postulated in simple potential 
scattering [13]. This postulate is necessary to prove that the S-matrix is uni- 
tary. Again, simply by inspection of Eqs. (9.4), we can see that the asymptotic 
states of the Rearrangement Model, according to these two formalisms, are 
not complete. Again, this point should not cause confusion. Our full states, 
Eqs. (6.4), (6.5), (6.6), and (6.7) are complete, as was shown in section VII. 



As a result, our "in" and "out" states form complete sets. However, the set of 
Haag-Ruelle type asymptotic states calculated according to either of the two 
formalisms is not complete. 

4. In potential scattering the Moller matrix, O, can be defined using the full in- 
teracting wave functions so that it is isometric even in the presence of bound 
states [1]. We see that the Haag-Ruelle asymptotic solutions, Eqs. (9.4), ob- 
tained by the use of 0, are certainly not orthonormal, whereas the original 
interacting wave functions were; therefore, the Moller matrix computed by 
their rules is not isometric, i.e. it is not norm preserving. However, the gen- 
eralized Moller matrix that we defined in Eq. (8.1) is not only isometric, but 
unitary. 

5. It is important to note that even though these asymptotic wave functions are 
neither orthonormal nor complete, they still lead to the correct S-matrix, as 
can be verified by calculating it using Eq. (9.6). If we had insisted upon the 
asymptotic wave functions being orthonormal and complete, we would have got 
the wrong S-matrix. 

6. Notice that because of this lack of orthonormality and completeness in the 
exact asymptotic states, the strong limits of Eqs. (2.26) are satisfied. Namely, 

lim (t) - V> in (<)] => 0, (13.1a) 
lim [*(<)- Vout 0> (13.16) 

lim Vin (0) ^ (0) EE fi(+Vin (0) , (13.1c) 

t^ — oo 

are automatically satisfied. This can be seen easily in the following way. For 
the first two equations above, the expression on the left hand side is just the 
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requisite full wave function, but with the delta function part removed. When 
we now take the norm and then take the limit, the remainder cancels giving 
zero. Similarly, the third equation above can be shown to be satisfied. 

On the other hand, if we had insisted that the asymptotic states are or- 
thonormal and complete, the wave function renormalization constants would 
have been missing from the asymptotic states. Thus, the delta function pieces 
would not have cancelled between the full and asymptotic states, and therefore, 
these pieces would contribute, and we would get a non-zero result. Thus, in 
this case, the strong limit would not hold. 

7. It is important to note that the reason that all these problems occur is that the 
full states are never considered. Most formalisms in the literature try to set up 
the problem in terms of the asymptotic states, and are thus forced to make as- 
sumptions regarding their properties and behavior. These assumptions are not 
necessarily correct in general, as is amply demonstrated by the Rearrangement 
Model, and other models such as the Cascade Model [7]. 

All this points out the importance of the correct normalization of the state vectors, a 
point already considered by DeWitt [31]. However, his work was restricted to the case of 
no bound states. The question of the correct description of the asymptotic states was also 
considered by Van Hove in his papers on the description of "persistent interactions" [4]. 
However, as noted in those papers, the formalism developed there does not deal with cases 
involving bound states, and does not deal with field theoretic scattering except for a few 
comments at the end. 
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In the multichannel case (such as rearrangement collisions), in the "channel Hamil- 
tonian" formalism, the statement is made that the basis states of one group of channels 
are not orthogonal to the others [10,11,13] because they are eigenstates of different free 
Hamiltonians. As we can see, in our model, the physical states C8<f), D<f>, and B6 are 
strictly orthogonal to each other. Evidently, this problem arises due to the use of "channel 
Hamiltonians" in the formalism. It is our belief that the method of splitting up the inter- 
action differently depending on which channel one is considering is fundamentally flawed 
because "every channel can be distinguished and is observable independently in experi- 
ments. This means that these channels should be orthogonal to each other" [32]. One 
method for ensuring orthonormality is given in [32]; however, this method still suffers from 
the flaws pointed out above. 

It is straightforward to see the problems caused by this lack of orthogonality. We are 
instructed, in these formalisms, to begin with asymptotic states. Let us first consider the 
"channel Hamiltonian" formalism. Then, the asymptotic states are the eigenstates of the 
"channel Hamiltonian" in the sector we are considering. As an example, let us consider 

\M B e(u))) -+ \M D <f>(u))), (13.2) 

where Mb is the physical B particle and Mjy is the physical D particle. We immediately 
notice, even before we consider any scattering, that the \MbO{lo))) state is not orthogonal 
to the \Mj)(f>(is))) state, as can be seen by inspection of Eqs. (5.10), (5.11), (5.12), (5.13), 
and (5.14). In other words, two experimentally distinct channels are not orthogonal to 
each other. This will clearly lead to the wrong S-matrix elements because it says that even 
if there is no scattering, there is a non-zero probability that the \MbO{lo))) state will turn 
into the \Mucf)(v))) state. We cannot even argue that the two states are "asymptotically 
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orthonormal" [10] because they clearly are not. This can easily be seen by observing that 
both \MbO{lo))) and \Mucf)(v))) have expansion coefficients in the "bare" \C6{u)(f>(v)) sec- 
tor. Therefore, as these states are neither orthonormal nor complete, we cannot have an 
isometric or unitary S-matrix, since orthonormality is necessary for isometry, and com- 
pleteness for unitarity. However, we have constructed a set of orthonormal (and complete) 
solutions for our system, a feat that many authors [33] tacitly assume is not possible, and 
have a perfectly isometric and unitary S-matrix. 

These problems with the S-matrix can be verified by explicit calculation. Since the 
calculation is tedious, we describe the method, and leave it to the interested reader to 
verify the results. Our interest is in the scattering of physical states, and so we must start 
by re-expressing the Hamiltonian, Eqs. (4.3) and (4.4), in terms of the operators which 
create the physical B and D particles. We denote these operators by B and Z>, respectively. 
They are found by inspection of Eqs. (5.10), (5.11), (5.12), (5.13), and (5.14), which are the 
wave functions for the physical particles. To find the expressions for these operators, we 
promote the states \C(f>{v)), \C8(lo)) } \B), and \D) to operators, all acting on the vacuum, 
and read off the expansions for the operators B and T>. In other words, 



We re-express the Hamiltonian in terms of these operators, which can be split into various 
channel Hamiltonians, from which the S-matrix is calculated. 

We can go even further than this. Consider the state |Mg#(u;))), which is a product 
state of the physical B particle and a free 6 particle. If we wanted the asymptotic state 




(13.3a) 



(13.36) 
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corresponding to this then, by the Haag-Ruelle protocol, we should find that the compo- 
nents of this state are only in the \C6cf>) and \B6) sectors, with no admixture of the \Dcf>) 
state. However, we can use our exact solutions to calculate this asymptotic state. We will 
find that this assertion will not hold true. 

To calculate the asymptotic state, we take the limit 

lim e lHct e- lHt \M B 6(uj))). (13.4) 

t^ — oo 

Inserting a complete set of states, we have 

/ dEdn lim e lHct e- im \E,n)){{E,n\M B 6(uj))) 
+ [dEdn lim e lHct e- im \E)) DD {{E\M B 6(uj))) 
+ [dEdn lim e lHct e- lHt \E)) BB {{E\M B e{to))). (13.5) 

J t^-oo 

Expanding each of the physical states, \E,n)), \E)) D , and \E)) B in terms of the bare states 
\C6{u)(f>(v)) , \B8(uj)) } and \D(f>(v)), and taking the limit, it is immediately obvious that 
the expansion coefficients in the \Dcf>(v)) sector are not zero. 

As a physical example, consider the case of a proton bound to a fixed nucleus by a 
potential Vp, and bombarded by a neutron which interacts with the proton and the nucleus 
through the potentials Vpjsf and V/v, respectively [32]. The total Hamiltonian of the system 
is 

H = Kp + K N + Vp + V N + V PN , (13.6) 

where Kp and Kjy are the kinetic energy of the proton and the neutron, respectively. The 
initial state, denoted by is given by 

= (13.7) 
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where 



H = H 1 + V 1} 



(13.8a) 



Hi = K P + K N + V P , 



(13.86) 



Vi = V PN + V N . 



(13.8c) 



Therefore, the initial state, is a product of a bound proton, cf>p(E l ), and of a free 

neutron (represented by a plane wave), up;(E l — Ef) } where Ef is the binding energy of 
the proton. 

Several possible reactions can occur giving rise to different final products. Let us 
consider four such reactions. 

1. Elastic or inelastic collisions. 

The proton remains bound to the nucleus, and the neutron is free after the 
collision. Therefore, the Hamiltonian is divided in the same manner as above. 

2. Exchange scattering. 

The neutron knocks out the bound proton and becomes bound to the nu- 
cleus. The Hamiltonian is then divided as: 



H 



H 2 + V 2} 



(13.9a) 



H 2 



K P + K N + V N , 



(13.96) 



V 2 



V PN + V P . 



(13.9c) 



Therefore, the final state is 



H 2 $ 2 j = -E/$2,/, 



(13.10a) 




(13.106) 
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3. Ionization. 

The neutron knocks out the bound proton and both are free after the 
collision. The Hamiltonian is then divided as: 

H = H 3 + Vz, (13.11a) 

H 3 =K P + K N , (13.116) 

V 3 = V PN + V P + V N . (13.11c) 

Therefore, the final state is 

H 3 $ 3 ,f = Ef$3,f, (13.12a) 

$ 3 J = up (Ef ) u N (E f - Ef ) . (13.126) 

4. Pickup. 

The proton and the neutron become bound and form a deuteron. The 
Hamiltonian is then divided as: 

H = H 4 + 14, (13.13a) 

# 4 = K P + K N + V PN , (13.136) 

V 4 = V P + V N . (13.13c) 



Therefore, the final state is 



# 4 $4J = £/$4J, (13.14a) 
$ 4J = u c (X, E f - Ef) <t> B PN (r, Ef) . (13.146) 



Here, X is the center of mass coordinate of the deuteron, and r is the internal 
coordinate of the deuteron. 
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Property 


Conventional 
Formalism 


Rearr. 
Model 


Asymptotic states normalized? 


Yes 


No 


Asymptotic states orthogonal? 


Yes 


Yes 


Asymptotic states complete? 


Yes 


No 


isometric? 


Yes 


No 


S-matrix unitary? 


Yes 


Yes 


Strong limit satisfied? 


No 


Yes 


He used? 


No 


Yes 


Additional property for the multiple-channel case 


Phys. states orthog.? || No 


No Yes 



Table 13.1: Comparison of the properties of the Rearrangement Model to 
various scattering formalisms. 

The final states given by Eqs. (13.10), (13.12), and (13.14), are eigenstates of different 
free Hamiltonians. Thus, in general, they are not orthogonal to each other, and the 
concomitant problems follow. 

The reason that these methods do not work properly is that the basis used is one 
in which bound-state eigenfunctions of the Hamiltonians that bind each fragment are 
multiplied by plane waves for the fragment motion [13]. In our model, because we have 
made no breakup, we get the physically reasonable result that the wave functions of the 
bound states are always orthogonal to the scattering states, and that the basis states 
of different channels are explicitly orthogonal to each other. We do not have to worry 
about making the explicit assumption that as the separation between the fragments goes 
to infinity, the overlap becomes negligible. This assumption may or may not be true, and 
leads to the problems with "persistent interactions" considered by Van Hove [4]. 
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We compare the results from the conventional formalism with those from the Rear- 
rangement Model in Table 13.1. 

In addition, even when it is not stated explicitly in the literature, it is often assumed 
that the spectra of the bound states and the scattering (continuum) states do not overlap. 
However, it is possible to construct models in which the spectra of one or more bound 
states overlap with the continuum [34,35]. Therefore, this assumption is not necessarily 
true, and will in general depend upon the details of the model under consideration. It 
is also possible to construct two different potentials which can lead to the same S-matrix 
with, in one case, redundant poles unnecessary for completeness, and in the other case, 
with the same poles being absolutely necessary for completeness [36,37]. This points out 
the need for resisting the temptation to identify the poles of the S-matrix with physical 
bound states of the system. 

More importantly, no authors have as yet worried about the evident normalization 
problem with the asymptotic states because they are always assumed to be normalized. 
These states are not normalized in the Rearrangement Model, and consequently, assum- 
ing orthonormality of the asymptotic states, in general, is very dangerous. In addition, 
we notice that in this model even though the asymptotic states are not normalized, the 
interacting states are. 

One approach that tries to avoid all these problems, especially in the cases of unstable 
particles and bound states, is that of analytic continuation [7,19,23,24,38,39] of the state 
space 7i into a generalized vector space Q. This has already been done for the case of the Lee 
model by Sudarshan, Chiu, and Gorini [19], Parravicini, Gorini, and Sudarshan [20], and 
by Bohm [40]. For instance, with this method, one can identify resonances and redundant 
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poles, and study the decay of a metastable quantum system. It can also be used for many 
other things, such as studying the Khalfin observation that the decay of a metastable 
system with an energy spectrum bounded from below can never be strictly exponential [41]. 
See the above references for details. 

14. Summary and conclusion 

In this work, we constructed a model that allows rearrangement collisions. We ex- 
plored the spectra and the complete set of orthonormal (ideal) eigenfunctions of this Re- 
arrangement Model in the Rearrangement Sector. Because of the structure of the effective 
Hamiltonian in this sector, we were able to solve the model exactly. In a similar fashion 
as for the Cascade Model [7], we find that the spectra can be interpreted as a B particle 
with energy < coupled to a 6 particle with energy lo } < to < oo; & D particle with 
energy Mp < coupled to a 6 particle with energy i/, < v < oo; and a C particle of 
energy coupled to 6 and <f> particles with energies u and i/, < u, v < oo. We see that 
the interacting field theory has a particle interpretation. 

Both the B and the D particles suffer mass renormalizations, and these mass renor- 
malizations alter the threshold of the B6 and Dcf> continua, respectively. In Eqs. (6.5b) 
and (6.6c), we also see the presence of both the mass and wave function renormalizations 
of the B and D particles in the plane wave parts of their respective wave functions. 

We have throughout emphasized the importance of using the comparison Hamiltonian 
(the diagonalized form of the effective Hamiltonian) because it is isospectral with the full 
Hamiltonian. Its spectrum differs from that of the free Hamiltonian by the alteration of 
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the B6 and Dcf> continua, and by the addition of a discrete A state. These effects are non- 
perturbative and, as emphasized in [7], can only be handled by a renormalized perturbation 
scheme in which Hq } not Hq } is taken as the starting point. 

Our results are surprising when compared to what we would expect from conventional 
scattering theory. We find that while the interacting state vectors are normalized, the 
asymptotic states are not. Moreover, the asymptotic states are neither orthonormal nor 
complete because of the presence of the wave function renormalization factors in the phys- 
ical Dcf> and B6 sectors. We note that this lack of orthonormality and completeness is 
absolutely necessary. If we construct the S-matrix from these states, we get the correct 
result (i.e. it is the same S-matrix as the one constructed from the full state). On the 
other hand, if we didn't allow the wave function renormalization factors because of our 
demand that the asymptotic states be orthonormal and complete, we would get the wrong 
S-matrix. Furthermore, for the strong limits Eqs. (2.26) to hold, we must again make sure 
to have these non-orthonormal and non-complete states. We also find that our physical 
C8<f), D<f>, and B6 states, while being the basis states for different channels, are strictly 
orthogonal to each other. Further, the Moller matrix, as defined in the literature is not 
isometric: it does not preserve the norm of the states. However, we defined a generalized 
Moller matrix which is not only isometric, but unitary. All these results are contrary to 
the usual formalisms of quantum scattering theory. 

More generally, we argued that the correct procedure, for any Hamiltonian, H , is to 
take its complete set of eigenstates, and an associated isospectral comparison Hamiltonian, 
Hq- The matrix of normalized eigenfunctions of H constitutes the generalized Moller 
matrix, which is unitary and intertwines H and Hq- 
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This model is a very simple one. However, even this simple model is enough to show 
the problems with conventional perturbation theory, and the conventional formulations of 
scattering theory. It is clearly necessary in the light of this model, and previous work on 
the existence of redundant poles in the scattering amplitude [36,37] and the presence of 
discrete solutions degenerate in energy with the scattering continuum [34,35], that a funda- 
mental re-examination be made of some of the postulates and assumptions of conventional 
quantum scattering theory. 
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A. Some Useful Formulae 

The following formulae are very useful for the calculations in the main text. By our 
definitions in section 6 we have the following ranges for our variables: 



< A < oo, 



(Al) 



< n < oo, 



(A2) 



< n < oo, 



(A3) 



with E being free to run over all values. 



We then have the easily proved identities 



</(A) 



2 



1 



T [/3(A)-/3*(A)], 



(A4) 



2vrz 
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\f(v)\ 2 

|g(A)l 2 

l/3(A)| 2 

IJ(^)I 2 

\a (//) | 2 

\g(E-\)\ 2 



1 

2vrz 
1 
2vrz 
1 

~2vrz 
1 



a (//) — a* (//)] 
1 1 



/3(A) /3*(A). 
1 1 
(/i) a*(/i) 



\P(E - A) | 2 2vrz 

i/(£-^)i 2 _ i 

|a(£-/i)| 2 2vrz 



1 



1 



P{E-\) p*(E-\) 
1 1 



Z B 6(E- \-M B ), 
■Z D 8(E-n-M D ) 



(A5) 
(A6) 

(A7) 

(AS) 

(A9) 



a (-E — /i) a* (i£ — /i) 
The equations (A. 8) and (A. 9) follow because E — A and E — fx can be less than zero, and 

thus pick up singularities at M B < and Mp < 0, respectively. On the other hand A and 

/j, are always greater than or equal to zero, and so cannot pick up any singularities. 

Another useful identity is 

1^7 (A) | 2 1 , Z B 



l(E) 



f ]\ |J v > 1 i_ 

|/3 (A) | 2 a (E- A) a(E-M B ) 



(AlOa) 
(A106) 



|a(/i)| 2 P(E - fj) P(E-M D )' 
We can easily show this by means of the contours in Figures 1 and 2. If we convert the 

integral in Eq. (A. 10a) into a contour integral by using Eq. (A. 6), we get 



dz 



+ 



Z B 



(All) 



2vrz/i Cl P{z)a{E-z) a (E — M B ) 
with the contour shown in Figure 1. Then, we make a change of variables from z to E — z 



to get 



~2vrz) ( 1} J Ci dz a(z)P(E-z) + a(E - M B \ 



(A12) 



with the contour shown in Figure 2. Now we deform the contour C4 and write it as the 
contour C3 plus the circle at infinity, while picking up the contributions from the residues 
of the integrand. Note that the circle at infinity gives no result, so we have 



■— )(-l)(-l)/ dz , s }, 

2m J K A J Jc 3 a(z)P(E-z) 
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+ (-^)«- 1) (- 1)(2 "»Mi^ + ^k)- (A13) 

The \[Z~b terms cancel, and the first two terms are Eq. (A. 10b), by definition. Therefore, 
Eq. (A. 10a) is equal to Eq. (A. 10b), and the identity is established. 
We can similarly show that 

1^7 (A) | 2 1 1 



dX 



\fl(X)\ 2 a(E - \)(\-u + ie) 

\f(p)\ 2 1 1 1 



a (/i) | 2 f3 (E - ji) (E - ji - v + ie) a(E - v) f3 (u) 



+ ^ ^ . (A14) 

P(E - M D )(E - M D - v + ie) a (E - M B ) (E - M B - to + ie) K } 

Using Eqs. (A. 7), (A. 8), and (A. 10) we can get another useful formula 

\f(n)\ 2 \g(E-n)\ 2 



dn 



a(n) \ 2 \f3(E - n) | 2 



_ 1 (E)- 1 *(E) |ff (E — Mj) | 2 \f(E-M B )? 

(-271-0 \/3(E - M D ) | 2 B \a{E-M B )\ 2 ' ' 

Finally, in Figure 3, we display the branch cuts and poles of which are used in 
showing the completeness of our solution set. 
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ure Captions 

Figure 1. Contour for Eq. (A. 10a) 
Figure 2. Contour for Eq. (A. 10b) 
Figure 3. Contour for the function 



